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PREFACE 


An  infinite^imensional  linear-programming  problem  pro¬ 
posed  by  W.  Prager  is  solved  in  this  Memorandum,  which  involves 
basic  mathematical  research.  The  problem  arose  in  an  elastico- 
plastic  structural-design  context.  is  hoped  that  the  solu¬ 

tion  of  this  problem  will  provide  insight  into  the  solution 
of  many  others  of  the  same  general  character. 


SUMMARY 


For 

subject 


the  problem  of  minimizing  the  Integral 
to  the  constraints 


-  f(x)  xg(y)  f{x)  If  0  X  _<  y  <  1, 

—  f(x)  xg(y)  —  X  +  y  f(x)  If  0  ^  y  ^  x  ^  1, 

the  author  exhibits  solutions  and  proves  both  that  they  satisfy 
the  constraints  and  that  they  have  the  extremlzlng  property. 

The  problem,  proposed  by  W.  Prager,  arose  In  an  elastlco- 
plastlc,  structural— design  context. 
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A  LINEAR  PROGRAM  OF  PRAGER'S 


1 ,  statement  of  the  problem  and  of  its  solution 

In  a  recent  letter  to  R.  Bellman,  V/.  Prager  proposed  the 
following  problem: 

It  Is  required  to  find  real-valued,  nonnegative  Integrable 
functions  f,  g,  on  [O,  l]  so  as  to  minimize 


subject  to  the  constraints 

-  f(x)  <  xg(y)  <  f(x) 

-  f(x)  <  xg(y)  -  X  +  y  <  fCx) 


lf0^x_<yili 
If  0  <  y  <  X  <  1. 


The  foregoing  problem  arose  In  a  structural-design  context 
of  elastlco-plastlc  type.  Although  the  writer  Is  not  familiar 
with  the  details  of  the  model.  It  Is  hoped  that  the  solution 
of  the  present  problem  will  provide  Insight  toward  the  solution 
of  an  entire  class  of  problems  of  the  same  general  character. 

The  solution  (f^,  g^)  given  below  Is  not  the  only  solution 
to  the  problem,  as  the  subsequent  analysis  Indicates.  It  turns 
out  that  gQ  can  be  perturbed  slightly  In  Its  tall  (beyond  7^/2) 
without  altering  feasibility.  Moreover,  since  the  other 
component  (fQ)  Is  defined  only  In  terms  of  the  values  of  g^ 
on  the  interval  [0,75/2],  Integrablllty  of  gQ  beyond  7^/2  is 
inessential.  It  Is  thus  seen  that  there  are  2°  solutions  of 
the  problem,  where  c  Is  the  cardinality  of  the  continuum. 
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Finally  we  remark  that,  corresponding  to  an  optimal  (f,  g) 
pair,  if  r  ^  f,  with  equality  holding  almost  everywhere,  then 
(r,  g)  is  an  optimal  feasible  solution  to  the  problem. 


Having  dispensed  with  this  digression,  we  now  contend 
that  the  following  (f^,  g^)  is  a  solution  to  Prager's  problem: 


In  Sec.  2  we  shall  show  that  (fQ,  g^)  satisfies  the  constraints 
imposed.  In  the  final  Sec.  3  we  shall  then  show  that  (f^,  g^) 
minimizes  the  objective  integral  over  all  pairs  (f,  g)  satis¬ 
fying  the  constraints. 

2.  VERIFICATION  OF  FEASIBILITY  OF  SOLUTION 

It  is  clear  that  the  given  constraints  can  be  rewritten 
as  follows: 

g(x)  >  0,  (A) 

f(x)  ^  0,  (B) 

f(x)  >  xg(y)  if  0  <  X  <  y,  (C) 

f{x)  >  xg(y)  -  X  +  y  if  0  y  <  X,  (D) 

f (x)  >  -  xg(y)  +  X  -  y  if  0  <  y  <  X.  (E) 
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We  remark  that  the  given  constraint 

—  f(x)  <  xg(y)  if  0  <  X  _<  y 

is  a  trivial  consequence  of  the  nonnegativity  conditionSj  (A) 
and  (B) . 

Before  we  show  that  the  given  pair  (fQ,  gQ)  satisfies  the 
constraints  (A)  ...  (e),  a  few  preliminary  remarks  are  in  order. 
Remark  1.  gQ  is  decreasing  and  convex. 

To  see  this  it  is  sufficient  to  show  that  on  the  set  where 
gQ  >  0,  its  derivative  is  negative  and  increasing.  But  on 
this  set,  we  have 


(1) 


go(x)  = 


/2  .  1 

/x  +  ^ 


It  is  immediate  that  the  right-hand  member  of  (l)  has  the 
requisite  properties. 

Remark  2 . 

A  numerical  calculation  gives 

go(0)  ^  -9407  <  1. 

Remark  3 »  Sq{^)  ^ 

This  is  an  immediate  consequence  of  the  two  preceding 
remarks . 


0. 


Remark  4 .  g. 


A  numerical  calculation  gives 


W)  •  ■ 


0593 


We  now  proceed  with  the  proof  of  feasibility  of  (fQ,  Sq) • 
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That  gQ  satisfies  the  constraint  (A)  is  obvious  from  the 
definition  of  gQ. 

Moreover,  it  follows  from  (A)  and  the  definition  of  fQ 
that  (b)  Is  satisfied  by  f^  on  the  interval  [O,  y2/2).  To 
see  that  Tq  is  nonnegative  on  the  remaining  interval,  it  suf¬ 
fices  to  show  that 


f^Cx)  =  -  XgQ 


if  X  > 


But,  at  X  =  y^/2,  the  right-hand  member  of  the  above  equation 
has  the  value 


—  Sq(0)  ~  ^ 

by  virtue  of  remark  2.  Consequently  (B)  will  be  established 
if  we  can  show  that  f  is  increasing  on  the  interval  [  y?/2,  l] . 
But,  since  the  function  —  1/2  maps  the  interval  (  y5'/2,  l) 
into  the  interval  (O,  ^/?/2),  and  satisfies  g^  >  0  on  this 
latter  interval  by  virtue  of  remarks  1  through  4,  we  may 
differentiate  fQ(x)  using 


Thus,  we  have 


3Gc 


But,  from  (l)  we  see  that 
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Substituting  this  value  In  the  above  formula  for  fQ  gives 

fQ(x)  =  1  -  go  -  ^  J  "  °  ^ 

by  virtue  of  remark  3.  This  establishes  the  desired  monoton¬ 
icity  of  fQ.  Hence  the  constraint  (B)  is  satisfied  by  (f^,  Sq) • 

To  see  that  (C)  is  satisfied  by  the  proposed  solution,  v;e 
must  prove  that 

fo(x)  >  xgQ(y)  if  0  X  <  y. 

But,  by  remark  1,  gQ  is  decreasing.  Hence  the  right-hand 
member  of  the  desired  relation  will  he  maximal  at  y  =  x,  and 
thus  it  is  sufficient  to  establish  that 

(2)  ^0^^^  -  ^  ^  * 

But  this  last  relation  holds  with  equality  on  the  interval 
[O,  /5/2),  by  the  definition  of  fQ  on  this  Interval.  Con¬ 
sequently,  we  need  to  show  only  that  the  relation  holds  on 
the  Interval  [ /5/2,  l] ,  i.e.,  from  the  definition  of  fQ  there, 
that 

fQ(x)  =  -  xgQ  -  !■)+■  X  -  yx^  -  ^  >  xgQ(x)  if '^  <  X  <  1. 

On  the  set  on  which  gQ(x)  =  0,  this  last  relation  holds  auto¬ 
matically  since  by  (b)  we  have  TqCx)  >  0.  Consequently,  we 
need  verify  the  above  only  for  points  for  which  gQ(x)  >  0 
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and  y2/2  x  _<  1.  But  at  x  =  y2/2,  we  have 

fo  ®o(°>  +4 


Upon  simplification,  using  v/ell— known  properties  of  the  loga¬ 
rithm  and  a  little  algebra,  we  obtain 

^0  ^  los  i/2  +  1)  . 


On  the  other  hand,  the  right-hcmd  member  of  the  desired 
inequality  yields,  at  x  =  /2/2, 


Thus,  the  desired  inequality  is  sharp  at  x  =  /2/2.  To 
establish  it  beyond  this  point  on  the  set  on  which  >  0, 
it  suffices  to  show  that  the  function 


fQ(x)  -  xgQ(x) 


-  xgQ 


XgQ(x) 


is  increasing  on  this  set. 

Upon  differentiating,  we  obtain 
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fo(x) 


XgQ(x) 


Xg 


0 


l^d  1 


4-  1  _  ^  /x2  _  1 

dx  2 


-  xgQ(x)  -  gn(x) 


By  the  previous  relation. 


I 

T 

=0 


-  1/2  /=  -  1/x, 


this  reduces  to 


f^(x)  -  xgQ(x) 


=  (l  -  j  “  xgQ(x) 


But  by  remark  3»  we  have  1  — 
have  —  xgp^(x)  >  0,  whence 


gQ(x)  >  0,  and  by  remark  1,  we 


fQ(x)  -  xgQ(x) 


>  0. 


Thus  (C)  is  established,  and  we  have  moreover  also  shovm 
that 


(2)  ^0^^^  -  ^SqCx)  for  all  x  e  [o,  l] 

To  establish  (D),  we  are  required  to  show  that 

fo(x)  >  xgQ(y)  -x+y  lfO<y<x  . 

By  remark  1,  Sq  is  convex.  Thus,  the  right-hand  member  of 
the  desired  inequality  is  convex  in  y  and  hence  is  maximal 
at  an  end  point  y  -  O  or  y  »  x.  Hence  it  suffices  to  establish 
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the  inequality  at  each  of  these  points.  But  at  y  =  0,  the 
right-hand  member  becomes  x  [^go(O)  -  1  <  0,  by  remark  2.  But 
since  TqCx)  >  0  we  see  that  the  inequality  is  satisfied  at 
the  lower  end  point.  At  y  =  x  the  inequality  reduces  to 

fo(x)  >  xg^Cx)  . 

But  this  last  is  the  inequality  (2)  already  established.  Thus 
condition  (D)  is  satisfied. 

Finally,  to  establish  (E)  we  are  required  to  show  that 

fo(x)  >  -  xg^Cy)  +x-y  ifO<y<x. 

By  remark  1,  is  convex;  hence  the  right-hand  member  of 
the  desired  inequality  is  concave  in  y.  Consequently,  if  there 
exists  an  interior  point  y^  e  (0, -x)  at  which  the  derivative 
of  the  right-hand  member  with  respect  to  y  vanishes,  this 
member  will  be  maximal  at  that  point.  Upon  differentiating 
with  respect  to  y  and  equating  to  zero,  we  obtain 

-  xgQ(y)  -1=0, 

or 


But  recall  that  this  equation  is  satisfied  by 

Vq  ~  ^  <  X, 

provided  of  course  that  /5/2  <  x  <  1.  Moreover,  y^  is  in  the 
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range  of  gQ  >  0.  Thus  for  x  e  [/2/2,  l]  it  is  enough  to  show 


On  this  interval,  however,  equality  holds  by  the  definition 
of  fQ.  Thus  it  remains  to  show  that  the  desired  inequality 
holds  on  the  Interval  0  ^  x  <  v/2/2;  i.e.,  we  are  required  to 
show  that 


xgQ(x)  >  -  xgQ(y)  +  X  -  y 


if  0  <  y  <  X 


Now  since  gQ  is  strictly  convex  on  the  set  on  which  gQ  >  0, 
and  M  -  1/2  is  imaginary  for  x  <  j2/2,  we  see  the  derivative 
of  the  right-hand  member  (with  respect  to  y)  cannot  vanish  at 
a  point  yQ  interior  to  [b,  /2/2).  Thus  the  right-hand  member 
is  maximal  at  an  end  point  y  =  0  or  y  =  x.  But  at  y  =  0,  this 
member  reduces  to  —  xgQ(O)  +  x,  and  we  are  required  to  show 


that 

xgo(x)  >  -  xgQ(O)  +  X 


if  0  <  X 


Since  this  inequality  holds  trivially  at  x  =  0,  w«  can  assume 
that  X  >  0,  and  this  last  desired  inequality  becomes 

^  ifO<x<N:^. 


But  since  gQ  is  decreasing  and  continuous  on  this  interval, 
it  suffices  to  establish  the  inequality  at  x  =y^/2. 

Prom  a  previous  calculation,  however,  we  have 
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'^"4}=  I  - 1  log  (y2  +  1), 

“  +  1  =  ^  ^  log  (^2  +  1)  . 

The  foregoing  inequality  is  thus  sharp  at  y  =  0.  Finally,  at 
the  upper  end  point  y  =  x,  the  Inequality 

Xgo(x)  >  -  XgQ(y)  +  X  -  y 


reduces  to  the  nonnegati; ity  condition  on  g^  alread 
Thus,  (E)  is  satisfied  and  the  verification  that  (f 
in  the  constraint  set  is  complete. 


y 

0’ 


established. 
Eq)  Is 


^ — PROOF  OF  optimality  of  solution 

The  proof  of  optimality  is  not  nearly  so  tedious  as  the 
rather  messy  analysis  of  Sec.  2,  which  showed  that  the  pro¬ 
posed  solution  (fg,  gg)  la  constraint  set.  The  reason 

for  this  unbalance  of  effort  Is  due  primarily  to  the  2-dllmen- 
slonal  character  of  the  constraining  Inequalities  and  the  fact 
(Which  we  do  not  prove  here)  that  equality  of  the  constraints 
can  hold  only  on  a  set  of  (2-dlmenslonal)  measure  zero.  This 
last  assertion  Is  an  aside,  however,  and  we  shall  turn  our 
attention  now  to  establishing  the  minimal  property  of  (f^,  g  ). 

L0  ^ 
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By  elemer.tary  calculus  and  a  numerical  computat  io: i ,  we  obtai 
J  =  J  ,  1  +  log(/2  +  1)  _y2  y-.ll68. 

Moreover,  we  have 

TqCx)  =  xgQ(x)  if  0  <  X  <  ^  , 


whence 


xgQ(x)dx. 


Also,  we  have 


Adding  equations  (3)  and  (4)  together  gives 


/.I 


fQ(x)dx  =  /  2xg  (x)dx  - 

{/r\  ^  iJ  rx 


^0  ^  l/Q 

Upon  making  the  change  of  variable 


I  /  p  T 

X  =  1 


in  the  second  integral  of  the  right-hand  member  of  the  above 
equation,  we  obtain 
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/  xg 

‘72  ° 


(x2  -  1]  dx 


/2 


=  ^  ^  x'SQ(x')dx'  . 


Thus,  the  integrals  Involving  cancel  and  we  are  left  with 
(5)  /'  f_(x)dx  =  J  . 


^0 


Next,  let  (f,  g)  be  any  pair  of  functions  in  the  constraint 
set.  Then  (C)  must  be  satisfied  in  particular  for  y  =  x  and 
0  <  X  <  j2/2.  Thus,  we  have 

(6)  f(x)  >  xg(x) 
whence  we  obtain 

(7) 


0  1  X  1  ^  . 


^0 


f(x)dx  >  /  xg(x)dx 


Moreover,  if  y5/2  <  x  <  1,  we  have  0  _  -  1/2  x,  so  that 

(E)  must  be  satisfied  if  y  =  /x^  -  1/2  and  /2/2  <  x  <  1;  that  is, 

(8)  f(x)  >  -  xg  -  i  j  +  X  -  /x^  -  i  if^^x<l. 

Integrating  this  last  inequality  over  the  interval  on  which 
it  holds  yields 


2  7'^ 

xg  \/x‘^  -  i  1  dx  +  J . 


/“I 

(9)  f(x)dx  >  xg  ^ 

Finally,  adding  (7)  and  (9)  gives 

f(x)dx  >  ^^xg(x)dx  xg  dx  +  j. 


But  upon  making  the  same  change  of  variable  as  before,  we  see 
that  the  integrals  involving  g  cancel  and  we  are  left  with 
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f (x)dx 


> 


J 


(from  (5) )  , 


v/hence  the  minimal  property  of  (f^,  is  established. 

A  final  closing  note  is  in  order.  In  Prager's  original 
statement  of  the  problem  it  v;as  tacitly  assumed  that  f  is 
Integrable  (\ia  the  objective  function).  No  such  restriction 
(as  we  have  made)  was  assorted  about  g.  Nevertheless,  a  proof 
of  optimality  with  this  restriction  relaxed  can  be  gi'.en  as 
follows : 

/  o 

If  v/e  replace  x  in  (8)  by  /x  +  1/2,  then  the  new  variable 
ranges  over  the  interval  [O,  /5/2]  and  we  obtain  in  place  of 
(8)  the  equl  alent  inequality 


Multiplying  this  last  Inequality  by  the  nonnegative  number 
x/./x"^  +  1/2  we  obtain 


Adding  the  inequalities  (6)  and  (ll),  which  now  hold  over  the 
same  range,  yields 


whence 
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Upon  making  the  appropriate  change  of  variable  in  the  second 
integral  in  the  left-hand  member  above,  and  the  same  change 
in  the  right-hand  member,  v/e  arrive  at  the  desired  inequality. 


In  accordance  with  an  earlier  assertion  about  the  multi¬ 
plicity  of  solutions  to  the  problem,  we  contend  that  the 
analysis  indicates  that  if  (f*,  g*)  is  any  pair  of  functions 
satisfying  the  following  conditions,  then  (f*,  g*)  is  a  solu¬ 
tion  : 

go(o)  _  g*(0)  i  1  +  go 

g*(x)  =  g^Cx)  if  0  ■  X  , 

max  ^0,  gQ  X  +  <  g*(x)  <  if  '^  <  x  <  1, 

f*(x)  ^  ^0^^^  with  equality  holding  for  almost  all 

X  e  [o,  l]  . 


(It  is  conjectured  that  these  are  the  only  solutions  to 
the  problem.  It  has  been  the  writer's  experience  in  such 
matters,  however,  that  a  proof  would  perhaps  be  rather  long 
and  will  not,  therefore,  be  attempted  in  this  paper.) 
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Graphs  of  these  solutions  showing  the  allowable  variation 
of  g*  are  appended.  Notice  that  there  is  a  discontinuity  in 
the  derivative  of  f*  at  x  =  . 
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